Polariton-polariton scattering in microcavities: A microscopic theory 
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We apply the fermion commutation technique for composite bosons to polariton-polariton scat- 
tering in semiconductor planar microcavities. Derivations are presented in a simple and physically 
transparent fashion. A procedure of orthogonolization of the initial and final two-exciton state wave- 
functions is used to calculate the effective scattering matrix elements and the scattering rates. We 
show how the bosonic stimulation of the scattering appears in this full fermionic approach whose 
equivalence to the bosonization method is thus demonstrated in the regime of low exciton density. 
We find an additional contribution to polariton-polariton scattering due to the exciton oscillator 
strength saturation, which we analyze as well. We present a theory of the polariton-polariton scat- 
tering with opposite spin orientations and show that this scattering process takes place mainly via 
dark excitonic states. Analytical estimations of the effective scattering amplitudes are given. 

PACS numbers: 71.36.-l-c,71.35.-y,71.35.Lk 
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I. INTRODUCTION 



The study of exciton-polaritons (polaritons) is a 
rapidly developing area in modern condensed matter 
physics. Cavity polaritons, which are two-dimensional 
excitons strongly coupled to trapped photons in quan- 
tum microcavities [J, exhibit a rich variety of nonlinear 
effects. These include stimulated scattering 2], polariza- 
tion rotation 0], bistable Q and multistable behaviors 
0, superfluidity [1, 0], and Bose condensation d, Q. 
Because of these unique properties, cavity polaritons are 
among the promising candidates for the implementation 
of future low threshold optoelectronic devices and their 
investigation thus has become topical. 

It is well established that the exciton-exciton interac- 
tion affects the exciton spin and momentum relaxation 
as well as the optical properties of semiconductors. The 
determination of the possible non-linearities governing 
optical properties is a fundamental question of prime im- 
portance in this field. Although it is quite clear that 
the excitonic interactions are the main source of non- 
linearity, the specifics of microcavity polaritons has not 
been entirely revealed and the polarization dependence 
of exciton-exciton has not been described in detail. The 
main complication of the theoretical description of the 
exciton-exciton interaction arises from the fact that ex- 
citons are composite particles, each being formed of two 
fermions: an electron and a hole. At first glance, such a 
problem seems quite similar to that of interacting hydro- 
gen atoms but the effective masses of positively and neg- 



atively charged carriers are comparable in semiconduc- 
tors; therefore a direct generalization of the techniques 
developed in atomic ph ysic s for the study of hydrogen- 
hydrogen collisions [13, Till . IT^ to the problem of inter- 
acting excitons is not possible and this forbids any treat- 
ment of the exciton-exciton interaction by use of an ef- 
fective Born-Oppenheimer (adiabatic) potential. More- 
over, even though the exciton-exciton interaction energy 
is much smaller than the typical exciton binding energy, 
it cannot be treated by standard perturbation theory be- 
cause of the electron-electron and hole-hole exchange pro- 
cesses. Indeed, it is impossible to know which electron is 
bound to which hole to form a given exciton and hence 
define the exciton-exciton interaction. This explains why 
much of the earlier attempts were based on oversimpli- 
fied models neglecting exchange interaction between exci- 
tons [1^ or the exchange of carriers between excitons 14 1 . 
This is not to mention the spin degree of freedom, which 
has been routinely neglected in the most part of previous 
works. 

Studies of excitonic nonlinearities are usually per- 
formed by using essentially different str ateg ies. One re- 
lies on the bosonization of excitons [Tsl . Il6l | an approx- 
imation that was believed to hold at very low excita- 
tion density. In this approach, excitons are assumed 
to be the only constituents of the dilute system under 
consideration and they also are considered as bosonic 
elementary excitations experiencing an effective repul- 
sion to avoid an overlap between the fermionic carriers 
wave functions. Under these assumptions, it was believed 
that by applying the Usui transformation (l7j the fun- 
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damental Hamiltonian could be mapped to an effective 
bosonic Hamiltonian. Despite unresolved issues at the 
heart of the bosonic approach (because of very restric- 
tive assumptions), its application has been justified by 
various groups [H, [H, [lo, HI], on the grounds of the 
numerical results they obtain, which describe some fea- 
tures of experimental data [23, lU, [2^, [2^, [S^ , and under 
the assumption that inclusion of a two-body exciton in- 
teraction term in the bosonic Hamiltonian is sufficient to 
account for fermion exchange effects. More recently, an 
effective spin-dependent exciton-exciton interaction po- 
tential of the Heitlcr-London type including long-range 
van der Waals terms was generated and studied [28j . 

Another strategy is based on the full fermionic treat- 
ment of the problem, which amounts to solve the second- 
quantized equations of motion [2^. Solutions of these 
equations rely on truncation schemes, which are tractable 
at the lowest order (Hartree-Fock) , i.e. assuming the 
density of the electron-hole system to be high enough 
to neglect excitonic correlations owing to the screening 
of the Coulomb interaction. However, at lower densi- 
ties, many-particle correlations have to be taken into 
account and the infinite hierarchy of equations satis- 
fied by multipoint correlation functions should be trun- 
cated at higher orders, which in practice proves quickly 
intractable. This technical difficulty originates in the 
choice of the strength of the Coulomb interaction between 
carriers as the "natural" physical quantity to character- 
ize and study the many-body problem in the interacting 
electron-hole system. In fact the Coulomb interaction 
cannot be treated as a perturbation when dealing with 
bound exciton states. Moreover, the random-phase ap- 
proximation, extensively applied to this type of problems, 
amounts to factorizc the multipoint correlation functions 
into two-point functions, neglecting in an uncontrollable 
way higher order Coulomb correlations such as biexci- 
tons. 

The optical excitation can also be considered as a per- 
tinent physical parameter for a microscopic theory ac- 
counting for all many-body correlations in the photoex- 
cited electron-hole system. Progress in the treatment of 
the equations of motion thus was made when a controlled 
truncation scheme based on powers of the optical excita- 
tion strength was developped [13, [U, [111 • For instance, 
in the contexts of four-wave mixing and pump-probe ex- 
periments in the low excitation regime, the optical non- 
linearities and related correlation effects are sufficiently 
well described at the third order of the applied laser 
field {x^^^ response). This method, called the dynamics- 
controlled truncation (DCT) scheme, allows to derive a 
set of closed equations of motion describing the dynam- 
ics of the exciton-exciton interaction at the mean field 
level and the four-body effects beyond mean field in a 
consistent way considering classical external laser fields 
[33I [11] . It also was succesfully applied to the study of 
six-wave mixing experiments with evidence of contribu- 
tions to the signals of x^^-' and x'^' processes [l^. The 
method was further improved with the inclusion of quan- 



tized electromagnetic fields in the Hamiltonian and the 
description of polaritonic effects [H, [13, [11] ■ 

The DCT scheme also provided a convenient frame- 
work for the analysis of biexcitonic correlations in terms 
of T-matrix |^]. More precisely, authors of Ref. (soj 
investigated the x*^'^'' optical response of a semiconductor 
quantum well to relate the many-body effects to the ex- 
citon scattering amplitude. However, in many cases such 
as, e.g., non-resonant pumping of the system, it is more 
convenient to work on the basis of excitons or polari- 
tons by introducing their effective scattering amplitudes. 
This approach proved extremely productive for the de- 
scription of polarization-dependent kinetic phenomena in 
microcavities [ll. Isl. l40j . 

The problem of the scattering of two composite ex- 
citons and related exciton-exciton interaction has also 
been addressed quite recently in the series of works by 
M. Combescot and co-workers lU, [H, [H [^, [H] , who 
developed a new technique based on commutation rules 
for the composite-exciton operators. Emphasizing the 
impossibility to define an interaction potential between 
two excitons because of the indistinguishability of their 
constituent carriers, their formalism was applied to study 
the scattering rates of excitons and certain differences 
with previous results were found [43]. The same group 
has considered the polariton-polariton scattering in mi- 
crocavities and proposed a new type of the optical non- 
linearity d^. 

The present work is focused on the detailed theoreti- 
cal description of polariton-polarion interactions in quan- 
tum microcavities. In Sec. II, we revisit the technique of 
Ref. j45| and present a somewhat simpler and physically 
more transparent description of exciton-exciton scatter- 
ing. We demonstrate how the bosonic stimulation of 
exciton-exciton (or polariton-polariton) scattering can be 
obtained in Sec. HI. We calculate effective scattering ma- 
trix elements for exciton-polaritons in microcavities and 
show how the previous results [Toj can be recovered. Spe- 
cial attention is paid to the nonlinearity caused by the 
saturation of the exciton oscillator strength in Sec. IV. 

Since an exciton-polariton possesses a spin degree of 
freedom, it is characterized by the projection of the total 
angular momentum of the electron-hole pair, -1-1 or — 1, 
on the growth axis. In Sec. V, we investigate in detail the 
effect of spin on the interaction of polaritons: it turns out 
that the scattering efficiency is strongly dependent on the 
mutual spin orientation of exciton-polaritons. We also 
present a microscopic derivation of the effective matrix 
elements for the scattering of polaritons with opposite 
spins, which is crucial for the linear polarization inversion 
observed in microcavities '4^. In Sec. VI, we present 
analytical estimations of the effective scattering rates. 



II. EXCITON-EXCITON SCATTERING 

In this section, we derive the basic quantities needed to 
calculate the scattering rates of composite bosons. To be 



3 



specific, we consider excitons in a direct band-gap semi- 
conductor quantum well. The splitting of the light-hole 
and heavy-hole bands is assumed to be large enough to 
neglect the population of the light-hole states. We first 
consider quantum well ground state excitons composed 
of electrons and heavy holes, ignoring their spin degree 
of freedom. The scattering rates derived here thus corre- 
spond to the scattering in the parallel-spin configuration. 
The scattering of excitons and exciton-polaritons with al- 
lowance for their spin is discussed in Sec. |Vl 

A. Wavef unctions 

Let (pi{r^,rii) be the wavefunction of a single exciton 
in the state z, where re, are the electron and hole two- 
dimensional position vectors. The state index « is a global 
index, which accounts for both the states of the relative 
motion of an electron and a hole in the exciton (such 
as Is, 2s, . . .) and the quantmn numbers of the center 
of mass motion, with center of mass wavevcctor K . We 
assume that the set of functions fi forms an orthonormal 
basis: 

J (p*{re,rh)(pi{re,rh)dredrh = Sij. 

In the second quantization framework, the exciton cre- 
ation operator i?| can be defined as 

Bj = ^,{K,k^,)albl, (1) 

where cpi{ke,kh) is the Fourier transform of (pi, and 

and 6^^ are the creation operators of the electron and 
hole with wavevectors and kh respectively. The se- 
quential action of the operators Bj on the vacuum state 
I vac) creates the corresponding numbers of excitons in 
the state i. For instance, the two-exciton wavefunction 

'if,j{re,,re,,rh,,rh,) = B|Bj|vac), (2) 

describes the pair of excitons in the states i and j. The 
wavefunctions obtained in this manner are properly an- 
tisymmetrized products of single exciton envelopes and 
can be considered as a zero-order approximation to the 
exact two-exciton states described by the Hamiltonian: 

Hexc ^Y.^a,+ 5]'K,v(r,, -nj, (3) 
ai ai,bi, 

where Tq, {a = e,h;l — 1,2) are the kinetic energy oper- 
ators of electrons and holes, and Va,b,, (a, b = e,h; I, V = 
1,2) the electron-electron, electron- hole and hole- hole 
Coulomb interactions energies. The notation in Eq. 



([3]) means that the terms with I = V corresponding to 
the same particles (a = h) should be excluded. 

In systems of interacting excitons, the basis of two (and 
more) excitonic states is overcomplete and hence non- 
orthogonal. Assume that there are N electron states and 
N hole states in the system. The total number of exci- 
tonic states is N^xc = N'^ (as one can take an electron 
and a hole in any of possible states). Now, let us calcu- 
late the total number of 4-particle states considering two 
electrons and two holes. For a pair of electrons (holes) 
one easily finds M = N{N — l)/2 states as two fermions 
cannot occupy the same state. Thus, the number of 4- 
particle states is ~ 7V^/4 for large N . However, 
the number of states with a pair of excitons defined by 
Eq. (O is iV2^c/2 ~ N'^/i since pairs (i,j) and (j, i) 
are equivalent. In a general case, all these pair states 
(i,j) have non-zero wavefunctions. Therefore, the num- 
ber of two-exciton states is approximately twice as larger 
as the number of 4-particle states. Thus, a two-fold ac- 
tion of the operator i?^ on the vacuum state generates a 
set of mutually dependent wavefunctions, and this yields 
two important consequences: (i) in order to calculate 
any transition rate one should orthogonolize the wave- 
functions; (ii) the total lifetime of an exciton due to the 
exciton-exciton interaction will not be equal to the sum 
of the transition rates of all possible two-exciton states. 

B. Effective matrix element of exciton-exciton 
scattering 

It is well established that excitonic nonlinearities 
in semiconductors are mainly caused by the exchange 
scattering of electron-hole complexes while the direct 
Coulomb interaction plays minor role. The exchange in- 
teraction has a short-range character (i.e. it becomes 
significant as the distance between the excitons' centers 
of mass becomes comparable to the excitonic Bohr ra- 
dius, as)- It means that e.g. for a two-dimensional sys- 
tem the exchange-induced corrections for a pair of ex- 
citons are governed by a small dimensionless parameter 
V = a?g/S 1 , where S is the surface area of the sample. 
If the exciton gas is described by its surface density n, 
the parameter v = na^. 

The overcompleteness of pair of excitons basis formed 
with the wavefunctions ^'(rei,'re2,r/i^,rft,2), defined in 
Eq. as basis vectors, calls for a special procedure 
to calculate exciton-exciton scattering rates. Let \f) 
be the initial and final two-exciton states, described by 
Eq. In general, these states are not orthogonal and 
we denote their scalar product as 

c = m- (4) 

Physically, it is obvious that the transitions should be 
considered only between the orthogonal states; otherwise 
part of the final state is admixed to the initial state. 
As we shall see below the constant C is proportional to 
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the parameter v and hence small. Therefore one may 
seek to consider orthogonolized pairs of states as slightly 
modified initial and final states, namely, 

\^')^\^)-a\f), |/') = |/) - (5) 

where a and (3 are some small complex coefficients 
(a,/? oc v). The scalar product of the states 
reads 

(z'|/'>=C- «*(/!/> -/3(zK). (6) 

The quadratic terms in a and (3 are neglected and the 
newly defined states are thus orthogonal provided that 
+ = C (the asterisk symbol denotes the 

complex conjugation). 

The matrix element of the total Hamiltonian between 
the orthogonalized states is 

(/'l-ffcxcIO = (/li^cxcN) -a(/|i7exc|/) -/3*(z|i/excN)- 

(7) 

Note that in the resonant situation where the ener- 
gies of the initial and final states are equal, the matrix 
elements (/|i?cxc|/)/(/|/) and (i|i/cxcN)/(«K> are equal. 
From Eq. ([7|) we obtain the transition matrix element 
Mt-,f as 

M.^f = (/'|i/cxcN') = (/|goxcN)-r ^''f?;^''\ (8) 

which is independent of the orthogonalization procedure, 
provided that the parameter v is small enough. Note, 
that in a general case the second term in Eq. ([8]) is not 
small as compared with the first one. Moreover, if one 
can formally separate the interaction part of the Hamil- 
tonian, Eq. ([5]) immediately gives the correct value for 
the transition matrix element as the matrix element of 
the interaction part. 

In order to establish the link between Eq. ([5]) and the 
results of Ref. [i^ we calculate the transition matrix el- 
ement in a simple case: two excitons initially occupy the 
same quantum state \i) = |00) while the final states for 
a pair of particles are defined by |/) = |12) = |21). It 



is assumed that the conservation of energy between this 
states is fulfilled. To proceed with the calculations the 
knowledge of three quantities is needed: the scalar prod- 
uct of the initial and final states, the matrix element of 
the full Hamiltonian between the initial and final states, 
and the mean value of the energy in the initial (or final) 
state. 

The scalar product of the states \i) and |/) given by: 

C = {f\i) = (12|00) = (vac I Si B2 I vac) 

--2a(?0),^ (9) 

can be derived using the commutation relations (45| or 
by expansion of Slater determinants. The quantity A ap- 
pearing in the equation above is the dimensionless Pauli 
parameter 44, 45] that describes an overlap between the 
initial and final state. More precisely, it is given by: 



"^[m i ) ^ J dreidre^dr^.dr,,, (10) 
V*m{re^ , ^/j^ )(/3; (r^^ , rh,)(p^ (Ve, , rh,)(pj (re, , J . 

Equation (O agrees with Eq. (7.4) of Ref. @|. It is 
worth noting that the integral in Eq. ^ describes the 
probability for two excitons to be at the same place and 
its magnitude is thus on the order of v. 

The matrix element of the Hamiltonian between the 
states |i) and |/) reads 

(/l-ffexcN) = (vac|SiS2i?excB^B^|vac) 

(u) 

where Eq is the energy of an exciton in the state 0. Here 
the quantities ^ and which respectively describe the 
Coulomb direct scattering, and the Coulomb exchange 
scattering between the in excitons, i.e. excitons in the 
states i and j in the terminology of Refs. 0, |4^, are 
defined as follows: 



^' TO i ) " / dT-eidre,dr,iidrft,¥3;,(rei,r,,J(p;(re2,rftJ^,(rei,rh.j95j(re2,r,iJ 



[Veeil 



^£2) + Vhhirhi - f/ij + Veh{ret -rh2)+ Veh{re 



(12) 



and 
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n J 
m i 



drei dre^ dr^^ dr,,, KX^ei , XC^^es , )ipt {r^, , r,,^ )ipj {r^^ , r^^ ) 

X [Vee{rei ~ r^^) + Vhh{rhi ~ T/ia) + yeh{rei - f/ia ) + ^^^(re^ - r^J] 

I 



(13) 



Finally, we need to calculate the mean value of the Hamil- 
tonian in the state We find: 



Combining Eqs. HT]) and (HH) we obtain 



1 



(14) 



(15) 



C. Transition rate 

Now we calculate the transition rate between the initial 
and final states. To this end, we use the Fermi golden 
rule and take into account that the states \i) and |/) are 
not normalized. The transition rate reads 



1 _2^|M._/| = 



ft m{f\f) 



5{2Eo-Ei-E2). 



(16) 



Since the matrix element M.i^f is proportional to the 
small quantity one can neglect the composite nature 
of excitons in the calculation of the normalization con- 
stants: (/I/) = (12|12) « 1 and = (00|00) « 2 [also 
used to derive Eq. p5|) ]. For the same reason one can 
use free exciton energies in the energy-conservation law. 
As a result we obtain: 



Here we show how to deduce the bosonic stimulation 
from the study of the exciton-exciton scattering matrix 
element. Following the approach developed in Ref. [ist . 
we consider the simplest case of a single exciton in the 
final state, thus the expected enhancement factor of the 
inverse transition rate 't~}^^ is equal to 2. The problem 
is complex because one needs to consider not two, but 
f/iree-particle states. 

The initial and final states are given by: 



\i) = |100) = B\{Blf\Y^c), I/) = |211) = B|(i?I)2|vac) 

(18) 

It is worth noting that, as in Sec. Ill Bi these states are 
neither orthogonal, nor normalized. As for the normal- 
ization constants we need only the main contribution: 



m « (/I/) 



(19) 



The calculation of the scalar products of the states 
|z) and I/) is more complex. As shown in Appendix iBl 
we obtain the scalar product {f\i) to the lowest (non- 
vanishing) order in v as: 



{f\i) = (vac|B^B2(i?J)'sl|vac) ~ -4A ( 2 q 



1 



(20) 



The matrix element of the total Hamiltonian between 
the initial and final states reads 



1 



li 



Si2Eo-Ei-E2) 



2 \ / 2 

1 J ^ \1 

' ' r-i ^^^^ 
which is equivalent to Eq. (11.18) of [4^ and Eq. (13) 

of Ref. [4^ , and we have thus proven that our simplified 
approach based on the orthogonalization of the initial 
and final states yields the same result for the exciton- 
exciton scattering rate as the more elaborate approach 
of Combescot and co-workers [491 ■ 



III. BOSONIC STIMULATION INDUCED BY 
EXCITON-EXCITON SCATTERING 

In the framework of the fermionic approach revisited 
here an important question arises as whether it is possi- 
ble to reproduce correctly the bosonic effects, mainly the 
bosonic stimulation of the exciton-exciton scattering. 



(/l^fcxcN) - (vac|BfB2i?cxc(ST)^S{|vac). (21) 

Rather tedious algebra (see Appendix [C| yields the lead- 
ing order contributions to Eq. pT|) : 



(/|i?cxcN) «rf+(Si+2So)(/|z), 



where 



2 
1 



2 
1 



(22) 



(23) 



is proportional to v. 

In agreement with Eq. ([5]) the transition matrix ele- 
ment reads 



M.^f = (/li/cxcN) - m'^^ « rf, (24) 
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and the exciton transition rate can be recast as 

(25) 

which is exactly twice larger than the transition rate with 
unoccupied final state, Eq. ^T7]i . 

Equation (|25|) clearly demonstrates that the fermionic 
approach yields, to the leading order in the bosonic 
enhancement of the scattering rate. It is worth noting 
that the corrections at the next order in i> reduce the 
transition rate as compared with Eq. ((25|) . In the case 
of N excitons in the final state, the procedure described 
above leads to the + 1 fold increase of the transition 
rate provided Niy <C 1. 

The transitions between fixed set of initial and final 
states thus can be described in the 'weakly-interacting' 
bosons approximation. We note, however, that such a 
procedure cannot be applied to the case where the tran- 
sitions cover all range of possible final states. Indeed, in 
order to calculate the exciton lifetime with respect to the 
transitions to all final states one needs to orthogonalize 
not only the initial and final states, but also all the final 
states between themselves. Since the basis of the pair of 
exciton states is overcomplete, it yields an extra factor 
of 2 in the transition rate, as shown in Ref. [i^ 



those of the electron and hole wavevectors fee, khi^^ clb^), 
yc cannot be neglected here in order to ensure the cou- 
pling of the photon and exciton with the same in-plane 
wavevectors. The light-matter interaction can thus be 
described by the following Hamiltonian 

i^cxc-ph - Y^^WPiA^ + W*P^A^^), (27) 

where is the photon creation operator and ly is a 
constant defined below. Note that the operators A^yt and 
A^ are true bosonic operators. 

In Ref. [i^, iJcxc-ph is defined with a light-matter cou- 
pling linear in the exciton operators: 

^fex"c-ph = E("^^-^- + (28) 

where il^ denotes the Rabi splitting, characterizing the 
linear coupling between a photonic mode and an exci- 
tonic state. Hereafter we consider only Is excitons and 
enumerate these states in accordance with their in-plane 
momentum. One can check that both Hamitonians, in 
Eqs. ((27)) and ((28)) . have the same single-particle eigen- 
states if the two constants W and fi^j satisfy 



IV. POLARITON-POLARITON SCATTERING 

Methods developed in Sec. II and III are applied here 
to the case of exciton-polaritons in semiconductor planar 
microcavities. Below we briefly review the basics of light- 
matter coupling in planar microcavities and afterwards 
calculate polariton-polariton effective scattering matrix 
elements due to (i) exciton-exciton scattering and (ii) the 
nonlinearity resulting from the saturation of the light- 
matter interaction strength. 

A. Polariton states in quantum microcavities 

Consider the interaction between electron-hole pairs 
with the light field. In the dipole approximation the for- 
mation of an exciton after the absorption of a photon 
is described by the dipole transition operator which can 
be conveniently expressed in terms of electron and hole 
creation operators as: 

= E AAu^Mu,^. (26) 

where h denotes the photon wavevector. The require- 
ment of momentum conservation in the process of photon 
absorption is satisfied by the presence of the Kronecker 8. 
Note that although the absolute value of the wavevector 
it involved in an optical transition is much smaller than 



nB.^W J ip*Jr,r)e-'"'-''dr. 

The cavity polaritons are formed as coherent superpo- 
sitions of the exciton states in a quantum well and photon 
states in a planar microcavity. In what follows we con- 
sider only polaritons arising due to strong coupling of the 
single cavity mode characterized by an in-plane wavevec- 
tor X with the Is heavy-hole exciton ground state with 
an in-plane wavevector K . The translational invariance 
in the plane of the structure imply the equality of the 
wavevectors of exciton and photon constituting a polari- 
ton, K = X. Moreover, we consider the states on the 
lower polariton branch only, therefore exciton-polaritons 
can be enumerated by a wavevector te and their creation 
operator can be written as 

Cl^a^Al + p^Bl (29) 

where unimportant extra quantum numbers are omitted 
for ease of notations. The prefactors and are Hop- 
field coefficients. 

The lower branch polariton dispersion reads: 

Eil = - (5 H 1 

-2V(^+2M^-2;^j 
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where S is the detuning between the photon and the exci- 
ton modes, Mx is the exciton total mass, and rrip^ is the 
photon effective mass, which depends on the width of the 
cavity and the refractive index. Here we have neglected 
the renormalization of the expression under the square 
root due to the broadening of the exciton and photon 
modes. Note that the Hopfield coefhcients can be di- 
rectly linked to the photon and exciton energies through: 



The matrix element of the total Hamiltonian H be- 
tween the states |z) and \ j) can be recast as a sum of two 
terms: one is related to the exciton-exciton interaction 
and is given, in accord with Eq. by 



2^ 



2 
1 



2C 



2 
1 



4<A 



2 
1 



2m. 



ph 



2M 



X 



(31) 



B. Effective matrix elements and 
polariton-polariton scattering rates 

Here we consider two main contributions to the 
polariton-polariton scattering cross section. One of them 
results from the exciton-exciton scattering, see Sec. |TT1 
the other one has its origins in the nonlinearity of the 
light-matter Hamiltonian, Eq. (P7|) . Below we calculate 
these two contributions, and we start by defining the two- 
polariton states as |x, x') = Cj^Cj^, |vac), which we write 
explicitely as: 



+ a^P^'A^Bl, + a^,l3^Al,Bl) |vac). (32) 



For simplicity of analysis, we consider the scattering 
of two polaritons initially in the same state \i) = \>Cot(o) 
to the final states |/) = \xiX2). Such a situation cor- 
responds to the optical parametric oscillators based on 
microcavities, which are widely studied both experimen- 
tally and theoretically, see e.g. Refs. (2|, |47|, El, Ull . In 
what follows we omit x in all places where it does not 
lead to confusion. The overlap of these states is 



The other contribution arises from the light-matter in- 
teraction part of the Hamiltonian and is related to the 
nonlinearity in the optical absorbtion caused by the sat- 
uration of exciton oscillator strength. It is also called 
^^photon-mediated exchange" (46| . The calculation of the 
matrix clement (/|i?exc-ph|*) yields the following result: 



(/|iJcxc-ph|i) = 

at/32*/3o' (vac| AiB2i?oxc-ph(st )2 |vac) + 
a;/3i*/32(vac|A2Si7fexc-ph(Si)2|vac) + 
2/3*/32*ao/3o(vac|BiB2i?cxc-ph45(l|vac) 



(35a) 
(35b) 
(35c) 



The physical sense of these processes can be illustrated, 
for instance, in the case of term Eq. (|35ap . The initial 
state of two polaritons |00) contains (with some weight) 
the product BgSQivac), but the pairs i?gi3g|vac) and 
vac) are not orthogonal and hence one of these ex- 
citons can recombine and emit a photon Ajjvac). After 
the emission, there is a probability for the remaining ex- 
citon to end up in the state 2. This state is nothing but 
a part of the polariton in CjCjIvac). Such a process is 
allowed by the momentum and energy conservation law, 
provided that 2Eq = E[ + and 2>to = Xi + H2. 
Using the identities given in appendix [Dl these matrix 
elements can be evaluated explicitly: 



{f\i) = /3r/32*/3o'(vac|i?ii?2(S^)'|vac) = -2(ilfi*f3l\ JJ 

(33) 

It is worth noting that the non-orthogonality of the po- 
lariton states is determined by the non-orthogonality of 
the exciton states because photons are true bosons. 

The total Hamiltonian of the system iJtot is the sum 
of three contributions: the excitonic Hamiltonian, i?oxc, 
given in Eq. ([S]), the light-matter coupling Hamiltonian, 
iJoxc-ph, given in Eq. (|27[) . and the photonic Hamilto- 
nian, 

ai ai.bii 

Y^iWPlA^ + W*P^Al) + J2 E^J'AlA^ (34) 



(/|i/exc-phN) = -2at/32*/3^P-^, (0, 0; 2)1^ (36) 

-2a^/3i*/32p_^^(0,0;l)H^ 

-2/3*/32*/3oao \p^, (1, 2; 0) + p^, (2, 1; 0)] W*. 

where 



p^{i,j;k) = J drdredrhip*{re,r)ip*{r,rh)ipk{re,rh)e"''^ 

(37) 

In the same fashion as above for the exciton-exciton 
scattering, we take into account the non-orthogonality 
of the two-polariton states Eq. ([5]) and use the modified 
Fermi golden rule Eq. (I16p , to obtain the polariton tran- 
sition rate as: 
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2 
1 



2 
1 



alP*Pt>P-><i (0- 0; 2)W - a*^Plf3tp-^, (0, 0; l)W 



-2f3lP*^(3oao [p-^o(l, 2; 0) + p_^„(2, 1; 0)] W* - PtP*2f3^o^ 
x5{2EP - - EP). 



2 
1 



7X 



E: 



f) \ao\^-aoP*o^R-a*f3on*i 



(38) 

2 



Equation ((38|) is the central result of this section. As 
compared with Ref. we use a different form of the 
light- matter coupling Hamiltonian, Eq. ^7} , which takes 
into account the fact that in the process of photon ab- 
sorption the electron-hole pair is created in the same 
point in real space. The analogous procedure with an 
effective Hamiltonian, Eq. (l28l) , would yield the same re- 
sult as published in Ref. [4^]. The main difference be- 
tween our results and those of Ref. [i^ is in the appear- 
ance of the p overlap integrals in our case, see Eq. ([57)1 . 
These overlap integrals contain products of three exci- 
tonic wavefunctions and one photonic wavefunction. This 
implicitly accounts for the fact that, e.g., in the process 
of photon absorption along with the creation of an ad- 
ditional exciton, the "resident" exciton in microcavity 
may change its quantum state. The non-linearity due 
to the light-matter interaction has been discussed also 
in Refs. [19l. [sol. [5l| . where overlap integrals of the form 
Eq. (p7|) also appear. However, the effects of the non- 
orthogonality of the exciton wavefunctions on the effec- 
tive matrix elements of polariton-polariton interaction 
were disregarded in these works. 

To conclude this section, it is worth noting that the 
contributions to the effective matrix element of polariton- 
polariton scattering due to light-matter interaction both 
in Ref. [i^ and in Eqs. ([55]) can be estimated as 

^ i^rin (see Eq. (IA6|) '). Thus, qualitatively, the differ- 
ence in the approaches of Ref. Refs. [l^,[50,[5l| and 
ours merely yields different numerical prefactors before 
the scattering amplitudes, while their dependence on the 
system parameters remains essentially the same. 

V. SCATTERING OF POLARITONS WITH 
OPPOSITE SPIN ORIENTATIONS 

So far, the effect of spin on the polariton-polariton 
scattering was not discussed and our consideration was 
restricted to a fully spin-polarized case: all polaritons 
were assumed to have the same projection of the momen- 
tum onto the growth axis. Below we extend our treat- 
ment to allow for the spin degree of freedom of polaritons. 

The spin projection of an electron on the microcav- 
ity growth axis is ±1/2 and the spin projection of a 
heavy hole is ±3/2. It is worth noting that the cir- 
cularly polarized cr"*" excitation creates an electron-hole 
pair (or exciton) in the spin state {e,hh) = (—1/2,3/2) 
which has a z-component of the total angular momen- 



tum TTiz = +1. The a~ polarized excitation creates an 
exciton in the spin state (1/2,-3/2) with z-component 
of the angular momentum being = — 1. The two re- 
maining exciton states with rriz = ±2, namely, (1/2, 3/2) 
and (—1/2, —3/2) are qualified as dark: they do not take 
part in the light-matter coupling but may serve as inter- 
mediate states for the scattering process. 



A. General remarks 

Here, we are interested in the situation where two ex- 
citons are created by two pump pulses with opposite cir- 
cular polarizations, cr"*" and . In such a case, a pair of 
excitons whose total 2-component of the angular momen- 
tum is equal to 0, is created. Let us introduce the exciton 
creation operator in the state with niz — +1, Bj as 

-B^+i = '^^(^<=''^'')«L,-i/2C3/2' (39) 

where (fi{ke,kfi) is the Fourier transform of the exci- 
ton envelope function (^(re,r/,,), ±1/2' ±3/2 ^-''^ 
the electron and hole creation operators. Here it is as- 
sumed that the subscript i denotes the orbital states of 
an exciton, say, i = {ls,K), and the second subscript 
denotes the z-component of the spin. Analogously, the 
exciton creation operator in the state with niz — — 1, 
B} _i, reads 

K-i= ^^(^^''^'')"L,+i/2^L,-3/2- (40) 

We need also two more exciton states which are dark, 
namely, those with = ±2. Their creation operators 
read 



bI+2^ Y1 '^*('^e,fe/0aI^_i/2&I^,3/2: (41a) 



^1-2= ^'('^'='^'')«L,-l/2^L,-3/2- (41b) 



9 



Let us now proceed with the pair of excitons states. 
Under the conditions described above the pairs of exci- 
tons are formed with opposite spin orientations: 

-1) = Bt^ii?[_Jvac), (42a) 



\i,+2;j,-2)=Bl^^Bl_,\vac) 



|z,-2;j,+2) -4_24+2|vac 



(42b) 



(42c) 



(42d) 



The following orthogonality relations greatly simplify our 
consideration: 

{i',m^;j', -m^\i,m^]j, -m^) = Su'Sjf , (43a) 



{i' ,rn^;j' ,-m^\i,-m^;j,'m^) =Su'5jj' . (43b) 

It means that within each set with a given |to^| = 1 or 2, 
the states are orthogonal. However, there is no orthogo- 
nality between the states with jm^j = 1 and \mz\ — 2: 



(^',+2;J',-2|^,+l;j,-l) 

= X! ^^"^ ' ^htWj' (fcei , fc/i2 (feei , fe/ii )<^j (fce^ , kh^ 



{k} 



(44) 



where, for ease of notations, {k} = k^^ , fega , k^i , kh^ . All 
other needed overlaps follow directly from Eq. 



B. Scattering via dark states 

We consider two channels for the scattering of excitons 
with opposite spins. The first scattering process appears 
in the Born approximation and is described by the matrix 
element 



M 



(1) 



(va.c\Bi,^+iB2',-iHBl^^Bl_^\vac). (45) 



We remind that the initial and final states \i) — 
bI ^iB\ _;^|vac) and |/) = B\, j^iB\, _-^|vac) are orthog- 
onal and hence an additional contribution due to non- 
orthogonality is absent in Eq. (|45[) . This matrix element 
can be calculated immediately, see Ref. (45j, Eqs. (8.25), 
(8.27) and (8.30): 



M 



(1) 



2' 2 
1' 1 



Equation (j46p is nothing but a direct Coulomb interac- 
tion between the excitons. In the case of cavity polaritons 
this contribution can be extremely small for polaritons as 
^(Q) cx [QubY , where Q is the transferred wavevector. 
In the case of cavity polaritons the matrix element of 
Eq. should be multipled by the product of Hopfield 
coefficients of initial and final states. 

Another possibility to have an interaction between ex- 
citons with opposite spins is to consider the scattering 
via intermediate dark states [s^l with \mz\ =2. In the 
case of excitons the splitting between the dark and bright 
states (determined mostly by the short-range electron- 
hole exchange) is usually small as compared with the 
typical energy of excitons; therefore such a scattering 
populates the dark states. This is not the case in mi- 
crocavities: the splitting between dark and bright states 
can be quite large (of the order of the Rabi splitting) and 
hence the dark states can play the role of intermediate, 
virtual states in the polariton-polariton scattering pro- 
cess. The related matrix element is weakly dependent on 
the exciton wavevectors, but it appears at the next order 
of perturbation theory. 

We carry out the calculation of the scattering rate be- 
tween the polaritons with opposite spins assuming that 
the splitting between polariton states and dark excitons 
states. A, exceeds by far the kinetic energy of polaritons, 
. This condition is fulfilled in quantum microcavi- 
ties under the resonant excitation of the lower polariton 
branch below or around the inflection point. 

The matrix elements of the interaction induced transi- 
tions from bright to dark states are needed. This process 
takes place due to the exchange of electrons (or holes) 
between excitons which results in the transfer of excitons 
from bright to the dark states. The Hamiltonian matrix 
element can be calculated using the fermionic approach 
as follows: 



(vac|Bi/,+2S2' -2-ffcxcSl,+i5| _;^|vac) = 
{El + E2){vac\Bi,^+2B2'.-2Bl+,Bl_,\va.c) 



nm 



^( m 1 ) (vac|Br.+2S2',-2S,l,+i<_i|vac) 

(47) 



(46) 



iE,+E2)x['^, n+r(;: I 



It is clear that such a process is solely determined by the 
exchange interaction. 

In order to calculate the matrix element of the exciton- 
exciton scattering via dark states we need: (i) to 
orthogonalize the intermediate states with respect to 
the initial and final states, and (ii) take into account 
the indistinguishability of the states |i, m^; j, — TOz) and 
|j, — m^; i, -l-mx). The first point can be easily achieved 
by the subtracting from the states of the pairs of excitons 
with \ mz\ = 2 the admixture of the states where \ mz\ = 1. 
The corresponding matrix element between the orthogo- 
nalized states differs from Eq. (|T7|) by the absence of the 
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first term (which is proportional to Ei + E2) [5g|. The 
second point can be dealt with by using the proper sym- 
metrization of the two-polariton wavefunction. In the 
case of opposite spin configurations one can introduce 
either symmetric (s) or antisymmetric (a) combinations 
as 



(48a) 



1 
71 



(48b) 

In the relevant wavevector range, the exciton-exciton 
scattering matrix elements weakly depend on the values 
of the wavevectors therefore the exciton-exciton interac- 
tion may be considered as short-range. In this case the 
interaction in the antisymmetric state is strongly sup- 
pressed as compared with the interaction in the sym- 
metric state because the overlap of the wavefunctions of 
excitons in the former state is small. 

In what follows let us consider the scattering of two 
excitons in a symmetric configuration, Eq. (I48ap . Let 
us denote \ii, 12] l)s the initial state, |/i, /2; l)s the final 
state. We denote as |ji, ^2; 2)s the states of pairs of dark 
excitons in a symmetric configuration which are orthog- 
onalized with respect to the initial and final states. The 
matrix element of the transition, calculated to the second 
order of perturbation theory, can be presented as (s^ 



Mflf = {E,,+E,,-E,,-E,,y^ (49) 



X (/i,/2;i|i7b-i,j2;2)(ji,j2;|i?Ki,i2;i). 

Here the summation is carried out over the different 
states (the states |ji,j2;2)s and |j2,ji;2)s are actually 
the same), and the subscript s is omitted. Using Eq. 
(HT)) this expression can be recast as 



Mflj = Y^{E,,+E,,~E,,~E,,)-' 



Jlj2 

out 



t I J2 «2 

ii h 



fi h 

m ( jl i2 
32 ii 



, (50) 



counterpart of Eq. (I50p for the antisymmetric wavefunc- 
tions differs by the signs in the square brackets only. We 
mention that since for free excitons in quantum wells, 
dark-bright splitting is usually very small as compared 
to the exciton kinetic energy, Eq. ((50|) cannot be used 
and the proper calculation of the matrix element to all 
orders of perturbation theory is needed. On the other 
hand, in microcavities the splitting between dark exci- 
tons and polaritons may be huge (of the order of the 
Rabi splitting) , which is why the extension of Eq. (|50p to 
the case of the cavity polaritons is perfectly applicable. 

This extension is straightforward: the energies for the 
initial states in the denominator should be polariton en- 
ergies, and energies of intermediate states are not af- 
fected by the light-matter coupling simply because these 
states are dark. Moreover the product of Hopfield factors 
/3n/3i2/3/i/3/2 should be introduced in Eq. Finally, 
the result for polaritons reads: 



M 



31 ,32 

t I f2 j2 \ I cin I fl 32 
fl 31 



t I 32 l2 

31 h 



f2 jl 

n I jl «2 
J2 ii 



(51) 



where the quantity is defined as 



Equation ([5T|l is the central result of this section. 
It is a microscopic derivation of the polariton-polariton 
scattering effective matrix element in the anti-parallel 
spin configuration, introduced in a number of works, see 
Refs. 0, Si- 
lt was demonstrated in Ref. [1^ that the negative rela- 
tive sign of matrix elements for polariton-polariton scat- 
tering in the parallel and anti-parallel spin configurations 
can lead to an inversion of the linear polarization in mi- 
crocavities in the non-linear regime. Provided that the 
main contribution to the polariton-polariton scattering is 
determined by the exciton-exciton interaction which, in 
turn, is determined by the repulsive exchange term ^™ 
(see Sec. IVip . the relative sign of the effective matrix el- 
ements is given by the sign of the energy denominator: 
{Ef^ + E^ - E^^ - E^). Obviously, if the dark states he 
above the bright ones (i.e. we consider the scattering of 
polaritons of the lower branch) the denominator is nega- 
tive which corresponds to an effective attraction of polari- 
tons and to the inversion of linear polarization. On the 
upper branch the dark states are lower in energy, there- 
fore the denominator is positive, the polariton-polariton 
interaction is repulsive and the inversion of linear polar- 
ization is not observed. 



n J 
m i 



3 n 
i ni 



VI. ESTIMATIONS OF SCATTERING MATRIX 
ELEMENTS 



and represent the Coulomb exchange scattering between 
the out excitons [3, EH]- It is worth noting that the 



Knowledge, even approximate, of the scattering ma- 
trix elements in contexts relevant to experiments is in- 
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structive. In this section we provide and discuss simple 
estimations obtained for both paraUel and opposite spin 
configurations. Detail of the evaluations of the integrals 
A and p^{i,j;k), respectively defined in Eqs. ^TU\i and 
([37)) . is given in Appendix A. 



A. Estimation of the scattering matrix elements in 
the parallel spin configuration 

1. Direct scattering 

As in Ref . [13] , we find that the direct term of the scat- 
tering matrix element, £,{Q), as a function of transferred 
wavevector Q, can be expressed as follows: 



K being the background dielectric constant. 
In the limit QaB <C 1, Eq. ([5^ reduces to 



(54) 



2. Exchange scattering 

Denoting ip{k) the Fourier transform of the 2D exciton 
ground state wavefunction Eq. (|A4p : 



(52) 



where 7e = me/Moxc, "fh = m/i/Mexc, with rUe, mh, 
Mexc = me + ruh being the electron, hole and exciton 
effective masses respectively; g{Q) = (1 + 0^(5^/4)^"^/^ 
is the Fourier transform of jV'lp)! , and Vq is the Fourier 
transform of the 2D Coulomb interaction potential: 



2aB 



S [l + (A:as)2]3/2' 



(55) 



we find that the Coulomb exchange scattering between 
in excitons, Eq. (jl3p . can be recast as 



(^?n _|_ 
^ee 



I (^in _|_ (^in _|_ /:in 



(56) 



Vq = 



kSQ ' 



(53) 



where 



k.q 



k,q 



2 

P+ Q 



r 



(57a) 
(57b) 



with 



P± ^ 2aiip + {ah ± ae)Q, 



and p is the wavevector characterizing the relative motion 
of excitons. 

The general expression of these four contributions to 
the scattering matrix element are quite complicated. 
However, taking into account that the wavevectors of the 
exciton-polariton center of mass motion are very small 
compared to the inverse Bohr radius, we can simply as- 
sume that the moduli of the vectors p, Q (and, therefore, 
P±) are negligible. It follows that 



Under the assumption that excitons are strictly two- 
dimensional so that Coulomb interaction is described by 
the potential Eq. ([55]) we obtain 



(59a) 
(59b) 







e2 

i 19 

Kas 


a% 
S 


'ieh " 


L f ™ ^ 


e2 

i 25 


a| 
S 



in agreement with Ref. [l9| . It is also instructive to com- 
pare the amplitudes of direct scattering and exchange 
scattering terms: 



a ^ai = ~J2^,v'(k + ^)^(k-^). (58b) 



k,q 



0.15(7^ -72)2(qas)3^ 1.5x10- 



(60) 



which shows that the exchange mechanism dominates the 
scattering process. To make the latter estimation we used 
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the following values: Q ^ 10^ cm"^ and Ub ~ 10~^ cm. 
Note that both direct and exchange scattering processes 
result in repulsive interactions of polaritons. 



B. Estimation of the scattering matrix elements in 
the opposite spin configuration 

The effective matrix element of the polariton-polariton 
interaction with opposite spins is the most difficult quan- 
tity to estimate. The problem stems from the fact that 
the dark states are involved and, as a consequence, the 
transferred momentum in the process of virtual transi- 
tion to the dark states can be arbitrary. The exact eval- 
uation of the effective scattering matrix elements will be 
a subject of a future work. Here we present an order-of- 
magnitude estmation of the considered effect. 

The initial and final states of an exciton are those with 
small momenta K,K' <C l/as . Therefore these mo- 
menta can be neglected and one can consider the scat- 
tering of two excitons from the ground state to the dark 
states with the momenta P and —P. Therefore, Eq. (|?T|) 
can be rewritten as 



\U{P)\' 



(61) 



Here hi^f is a product of Hopfield factors for initial and 
final states, cf. Eq. ([5T|l . A is the splitting between dark 
and bright states, and 



U{P) 



P 
P 



P 
P 



The summation over P can be transformed into an inte- 
gral as 



S Mx 



d£ 



(62) 



To make an estimation we assume the following model 
for U{e): U{e) = Uq for e < Eq and otherwise. An 
integration in Eq. ((62)) yields 



M 



(2) 



■Un In 



27r 2^2 



(63) 



We set Uo = e'^/{K,aB){a%/S) (see Eqs. ((5ga|) and (|59bp ). 
and Eq = / {Mxa\) (because the wavevector cut-off is 
at P ~ 1/as and, hence, the cut-off energy is h^K"^ /Mx)- 
Thus, one obtains 



M 



(2) 



"2 ^ ' a| Mxal 



In 1 



The numerical factor is omitted here. 

Depending on the relation between 
and A there are two limits: 



Mxaj^A 

(64) 

nyiMxal) 



1. if e 



2Mxa 



^ ^ A, we obtain 



M 



(2) 



1 



9 \ 2 9 



2. if e 



2Mxa 



A \ naB 
^ ^ A, we obtain 



s ■ 



(65) 



M 



(2) 



a 



\KaB J S fi2 



B Mxal 



Mxa%A 



(66) 



So far we have assumed that the light-matter coupling 
is relatively strong, so that the dark-bright states split- 
ting A, being equal to half the Rabi splitting, exceeds 
by far the matrix elements of polariton-polariton inter- 
action. In this case we can limit the calculation of the 
scattering amplitude in the opposite spins configuration 
to the second order of perturbation theory. The situation 
is different if the light-matter coupling strength is small. 

In order to gain a qualitative understanding of the sit- 
uation with small values of A, one can neglect the light- 
matter interaction entirely and consider the interaction of 
bare excitons. The exciton-exciton scattering via excited 
states is discussed in Ref. [l^l, where it is demonstrated 
that an interaction of polaritons with opposite spins can 
result from the scattering via such states as 2s, Ip, . . . An 
estimation given in [2CH has the same dimensional param- 
eters as Eq. 

Moreover, the condition Kgb ^ 1 shows that the 
scattering takes place in the so-called low energy region. 
Thus, the scattering rates can be strongly enhanced. Un- 
der the assumption that the Rabi splitting is smaller than 
the bi-exciton (bi-polariton) binding energy, one may use 
the results of a gene ral theory of two-dimensional scat- 
tering and write [SOllssl Is^ . [ssj 



MxS 



In 



e 



(67) 



Here e is the energy of the relative motion of a pair of 
excitons (or polaritons), i?o > is the bi-exciton (bi- 
polariton) binding energy. Note that scattering via bi- 
exciton or bi-polariton states [s^ may be even more effi- 
cient than the scattering via dark states. 



VII. CONCLUSION 

The scattering of excitons and exciton-polaritons in 
two-dimensional semiconductor microcavities has been 
considered here in the framework of the fermionic com- 
mutation technique Ref. [44]. Our results can be sum- 
marized as follows: 

1. The basis of two-excitonic states is overcomplete 
which requires an orthogonalization of initial and 
final wavefunctions for the calculation of the matrix 
elements and scattering rates. 
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2. For a fixed finite number of inital and final states 
the scattering rates calculated within the fermionic 
approach and bozonization scheme are equivalent. 

3. In the low-density regime the bosonic enhancement 
of exciton-exciton scattering is reproduced in the 
fermionic approach. 

4. An additional contribution to polariton-polariton 
scattering rate which arises from the non-linearity 
in light-matter interaction is derived. 

5. The scattering matrix elements for exciton- 
polaritons with anti-parallel spin configurations are 
calculated. 
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APPENDIX A: EVALUATION OF p AND A 
OVERLAP INTEGRALS 

To evaluate the integral p^(i^j]k) as defined in 
Eq. (|37p . we assume that the exciton wavefunction takes 
the following form: 



(^Jr^r/O^^e'^-^V-Cp), (Al) 

where Ka is the exciton center of mass wavevector, R — 
leXe + ^hfh is the position of the center of mass and 
p = Te — is the electron-hole relative coordinate. Here 
S is the normalization area, ■(/'(p) is the relative motion 
wavefunction. Substituting pi = — r and pi = r — r}^ 
and integrating over r we obtain 



dpi y"dp2-i/'(Pi)V'(P2)'0(pi + P2)x 



exp [ipi • (~7e>*i - Te-'^j + Te-J^o)] exp [-ip2 • [-^ih^i - IhXj + JhXo)]- (A2) 



r 



The characteristic length of the ijj function variation and introducing the dimensionless quantity 
is the excitonic Bohr radius, while typical values of 
the wavenumbers >fij-,o ^ 0,3^ because the polariton 
wavevectors are determined by the photon wavevectors. 
Therefore, in the integrations over pi and p2 one may 
omit the exponential terms, which yields: 



Vs 



(A3) 



X y^pi y dP2^/'(pi)V'(P2)'0(pi + P2). 



Assuming the following form for the wavefunction of 
the relative electron-hole motion in the exciton 



e ''^e ''^ exp 



we obtain 



Pi+ P2 + 2pi/52 cos 6 



P>t{hj;k) 



SW 



0.896, (A5) 



(A6) 



cxp(-p/aB) ; 



(A4) 



where W is defined in Eq. (|77|). 

In the same manner the A overlap integral 



A = y drei drea drh^ dr^, ip*^ {ve, , rh^)^*^ (re^ , r,, J(p, {re, , rhjfj {ve^ , J . 



(A7) 



reduces to the following expression 
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is 



J dpi J dp2 j dri;{pi+r)ijj{pi)i}j{p2-r)ip{p2), 



(A8) 



which can be further simpHfied: 



8ai 



^>Ci-\->tj.^j-rL-\-^rL*-^ ' 



(A9) 



where is given by 



J = J dpidp2 J dr j dOi J de2PiP2re-P'e-P' 

X exp —\Jp\+r'^ — 2/9ir cos 9i (AlO) 
X exp T o ^ 



P2+ r"^ + 2/92?' cos 62 



and is numerically evaluated: ~ 3.95. 



APPENDIX B: CALCULATION OF THE 
SCALAR PRODUCT INVOLVING THREE 
PARTICLE STATES 



Using the commutation rules outlined in Eqs. (6.32)- 
(6.35) of [4^, we can exchange the order of B2 and (i?o)^ 
in Eq. I|20p and, afterwards, swap B2 and B\. Finally, 
there are two types of terms that do not vanish: 



^2(vac|i?2BtBt|vac) 



n 
2 1 



+ A 



2 
n 1 



and 



^2(vac|B2BtB||vac)A 



71 

2 



(Bl) 



(B2) 



There is a significant difference between the magnitudes 
of the terms given in Eqs. (|Bip and (|B2p . The contri- 
bution (|Bip is small as compared to (IB2p because the 
latter is determined by the term with n = 1, where 
(vac bIb\b\\ vac) « 2, while the magnitudes of all the 
terms in (IBll) are further reduced because of the factor 



{-vac\B\B\Bli\vac) oc v. Retaining only term with n = 1 
in Eq. (|B2p we obtain the right hand side of Eq. (^0)) . 



APPENDIX C: CALCULATION OF THE 
MATRIX ELEMENT OF HAMILTONIAN 
INVOLVING THREE-PARTICLE STATES 

Performing the commutation of the operators i?cxc and 
{Blf in Eq. ini), we obtain four terms 

(vac I Bl B2 H,^,{bI ) 2 B 1 1 vac) 

{v£,c\BfB2{Bl)^H,y,,Bl\va.c) (Cla) 

+ 2Eo{Y(ic\BfB2iBl)'^Bl\v'Ac) (Clb) 

+ 2{vac\BfB2BlV^ Bl\vac) (Clc) 



EC 



n 
m 



(vac|i32i?2i?t^Bti?t|vac). (Cld) 



The first two terms on the right hand side of the equation 
above, which we denote a and b can be evaluated easily: 



2Eo{f\t). 



(C2a) 
(C2b) 



Term Eq. (jClcP can be recast in a similar to the term 
Eq. (jCldp form: 



2Ve 



n 
m 



Yac\BfB2BlBlBl\Yac) 



nm 

and the term Eq. (jCldp can be recast as 

d= (vac|B„S„i(s|)2BiB||vac)* 

+ 2{vac\BnB^Bl{l~ Dii)Bl\Ya.c)* 

k 1 
1 1 



2 2 A* 



vac|B„B™B]^B||vac) 



(C2c) 



(C3a) 
(C3b) 

(C3c) 



where Dn is the deviation- from-boson operator intro- 
duced in Ref. |44| . The leading order contribution (cx 
is given by the first summand in Eq. (|C3bp . It yields: 



d = 4e 



2 
1 



-4f 



2 
1 



(C4) 



One can check that Eq. (jC3ap gives a zero contribution, 
and Eqs. (|C3cP and (jC2cp yield similar expressions but 
with extra A factors, which give corrections of the order 
of ly^ to the effective scattering rate. 



APPENDIX D: MATRIX ELEMENTS OF 
EXCITON-PHOTON INTERACTION 
HAMILTONIAN 

To calculate the contribution of absorption nonlinear- 
ity to the polariton-polariton scattering rates, knowledge 
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of the overlap between the polariton state P^Bl\yac) and the two-exciton state -Bj-Slkac), is useful: 



(vac|BiB2P^st|vac) ^ 

/ dreidr^i / dre^dr^^^ / dreodr^j, 

fef ,fej^ fe2)'^2 '^O'^^O ^^'^'^ 

exp [ifefrei +1^1 ''ea +''^k2rh2 - iklVeo - ifcorho](/3t(rei,rhjv'2(»'e2,»'h2)</'o(r'eo.»'ho)x 

( vac I Ofee ttfee al,^ a^e I vac) ( vac 1 64,5. 6^^ 6]^^ | vac) = 

j dLrdre6.rh^pl{r,r)e^'*'^ip\{re,rh)^o{re,rh) + j 6.rdredrhV>\{r,r)e^'''^ipl{re,rh)v>o{re,rh) 

- j drdredr/,(p*(re,r)(p2(»',T?i)e'*"'v?o(re,r,i) - j drAre6.rhV>l{re,r)ip\{r,rh)e^'''^ipo{re,rh), 
The above result was found using the following identity: 

(vac|afeeafe|ai,^a^e|vac)(vac|6fc5.6fej6l,^6^Jvac) = {5ki,ki5^,,k-- 5k^^,k%5^„kl) x {^k'sM,^>'hM ^ ^KM^>'hM 



(Dl) 
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